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ABSTRACT 

We have investigated the cosmological implications of the bulk viscous cosmological model 
with variable G and A. These results are found to be compatible with the present obser- 
vations. The classical cosmological tests for this model encompass the Freese et al. ones. 
The model has some spirits of the Standard Model. The inflationary solution which solves 
the Standard Model problems is obtained as a special solution. The influence of viscosity 
is shown to affect the past and the future of the Universe. 
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1. Introduction 



We have recently discussed a bulk viscous model with variable G and A in a scenario 
which conserves energy and momentum [1]. We have found that the model has many 
interesting features. Various models in the present literature are shown to be equivalent 
to bulk viscous model. The bulk viscosity coefficient (77) is usually taken to depend on 
the energy density (p) of the cosmic fluid, viz. 77 = i] p n , where 770 is a constant and n is 
the viscosity 'index'. It is shown that only models with n = 1/2 are structurally stable. 
Beesham has shown that the Berman model is equivalent to a bulk viscous model. We, 
however, have shown that this equivalence is quite general. 

The introduction of bulk viscosity is usually accompanied with particles creation, since 
this term appears as an energy source in the energy conservation equation. It is also shown 
that the introduction of a variable cosmological constant solves the entropy problem. In 
this model neither of the two cases appear. This can be understood as the fact that 
the conservation of energy in cosmology is very restrictive. The conspiracy between the 
gravitational constant (G), cosmological constant (A) and the bulk viscosity is required 
to accomplish this. These three parameters are found to be related by a new equation. 

This model shares some spirits of the Standard Model on one hand and the inflationary 
model on the other hand. But, in this case the inflation is caused by viscosity effects. The 
notable feature of this model is the fact that the Standard and the inflationary Models are 
augmented in one model, and that the inflationary solution does not require the Universe 
to have the critical energy density. 

It is interesting to note that the inflationary solution can be obtained with or without 
the cosmological constant. Also note that this solution in this model requires 77 = rjop 
and H oc tjq 1 , where H is the Hubble's constant. So the smaller is the value of 770 the 
bigger is the inflation rate. This solution has been obtained by Murphy (1973) where he 
attributed the viscosity effect to gravitons production in the graviton-graviton scattering. 
The physical existence of this bulk viscosity is still needed to be studied. On the other 
hand, if it has been proved that the gravitational constant is varying with cosmic time, 
then this cosmology becomes the next candidate. 

A variation of G has many interesting consequences both in geology and astrophysics. 
Canuto and Narlikar [9] have shown that the G varying cosmology is consistent with 
whatsoever cosmological observations presently available. 

In this paper we would like to discuss the cosmological implications of this model. 
Recently, Waga [2] has discussed these implications in decaying vacuum energy flat models 
but with G constant. We have shown that the model of Freese et al. can be reproduced 
from this model. As for Freese et al, the nucleosynthesis constraints restrict (3 to be 
P < 0.1. Complying with the recent values of the age of the Universe our model fits very 
well. 
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2. The Model and Field equations 

In a Robertson Walker universe 

dr 2 = dt 2 -R 2 (t)[ - , - + r 2 (d9 2 + sin 9 2 d<j) 2 )] (1) 
1 — kr z 

where k is the curvature index. 

Einstein's field equations with time dependent cosmological and gravitational "constants" 

iV - \g^R = SttGT^ + A 9lu , (2) 

and the perfect fluid energy-momentum tensor 

T,!,, = (p + p)u^u v - pg^ (3) 
, in the presence of bulk viscosity, give [1] 

3(p* + p)R = -(^p + p + ^)R. (4) 

where p* = p — 3r]H and r\ is the coefficient of viscosity. The conservation of energy 
requires 

p + 3H(p + p) = . (5) 
The pressure (p) and the energy density (p) are related by the equation of the state 

p = (7 — l)p , 7 = const. . (6) 

Therefore, eq.(4) yields 

9 v HR=(^p+^)R (7) 

We take the ansatz [6,8] 

A = 3(3H 2 , 13 const.. (8) 

Lima et al. argued in favor of a term of the type A ~ H 2 . They assumed, on dimensional 
grounds, that A ~ J- (7 s -)™ where tn is the Hubble's time, tpi and £pi are respectively the 

t pl tu 

Planck's time and length. They take n = 2 in order to get rid of the % dependence, since 
the gravitational effects today are described by some classical theory of gravity. Therefore, 
since tu ~ H' 1 they concluded that A oc H 2 . On the other hand, the constancy of the 
parameter (5 = p v /(p + p v ) (p v is the vacuum energy density) proposed by Freese et 
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al. is equivalent, at the level of the Einstein's equations, to considering A = 3[3H 2 too. 
Equation (7) can be written as 

H G A' 

where a prime denotes derivative w.r.t scale factor R while a dot is the derivative w.r.t 
to cosmic time t. In what follows we will consider a flat universe, k = 
The density parameter is given by 

n = |g! = (i-/j). (io) 

The solution of eq.(9), using eqs.(5), (6), (8) and (10) yields 

R(t) = [3£> 7 (1 - n )][V37(l-n)] t [l/37(l-")] ? £ = const ( U ) 

The Hubble parameter is H = 3( - 1 i TO ) ^ and the deceleration parameter is g = lo — 
2 — 3n , < n < 1 ( the subscript denotes a present day quantity). This gives 
-1 < go < 2. 

3. Lookback Time 

The radiation travel time (or lookback time) t — t for a photon emitted by a source at 
instant t and received at to is 

f R °dR 

where i? is the present scale factor of the Universe. Using eq.(ll) one can write 

R = A't 1 / 3 ^ A' = const. , (13) 

it follows that 

f = l + z = ( t j)^\ (14) 

or 

t=(l + z)-^- n h . (15) 
Using eqs.(ll) and (15) we can write 

or 

Ho(t -t) = ^r^y[l " (1 + ^)- 3(1 - n) ] • (17) 
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For small z one obtains 



H (t -t) = z-(2-^n)z 2 + ... (18) 



and with q = 2 — 3n this transform into 

H (t -t) =z- (l + q/2)z 2 + ... , (19) 

If we take z — > oo in eq.(17) we readily obtain H Q t = g^z^ ; n — 1/2 gives the well-known 
Einstein-de Sitter (ES) result, 

H (t - t) = ?[1 - (1 + 3)-3/2] . (20) 

We recall that the FRW flat models, as preferred by inflation, have q — 1/2 and = 
2q — 1, whereas the dynamical estimates yield 0.1 < Q < 0.4 [5]. Such a problem can be 
resolved in our scenario. 

4. Proper Distance d(z) 

A photon emitted by a source with coordinate r = r\ at time t — t\ and received at time 
t by an observer located at r = will follow a null geodesic with (#, 0)=const. . The 
proper distance between the source and the observer is given by 

r R o dR 

i=R °L rr ' < 21 > 

hence 

rf = r lJ R = ^-[l-(l + z) 3 ™- 2 ] . (23) 

For small z this becomes 

H d = z- H(i _ n )z 2 + ... = z - 1(1 + g)^ 2 + ... (24) 

Equation (23) gives the Freese e£ a/, result for the proper distance if we set Q = 1 — (3 = 
2 — 2n viz., 

^ d=^A_[l-(l + ^-(3/^]. (25) 

As remarked before the ES result is a special case of eq.(23). Equation (23) shows that d 
is maximum for n — > 1 (de Sitter) and minimum for n — ■> 1/2 (ES). Now if n < 2/3, at to> 
the observer cannot see any source beyond the proper distance of the particle horizon [2] 

* = = & (26 > 

and for 2/3 < n < 1, d(z = oo) is always infinite. 
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5. Luminosity Distance 

This generalizes the inverse-square law of the brightness in the static Euclidean space to 
an expanding curved space [2] 

d L = (^) 1/2 = ri Ro(l + z) (27) 

where L is the total energy emitted by the source per unit time, £ is the apparent lumi- 
nosity of the object and r\ is the coordinate distance. For flat universe d = riR (l + z) 
and therefore one can write 

d L = d{l + z). (28) 

Using eq.(23) one gets 

H d L = ±±^[1 - (1 + z)^} . (29) 
For n = 1/2 the above equation gives back the ES result, 

d L = ^[(l + z)-(l + z)^]. (30) 

For small z eq.(29) gives 

H d L = z + l(l- q)z 2 = z + i(3n - l)z 2 + ... (31) 

6. Angular Diameter 

The angular diameter of a light source of proper distance D at r = r\ and t — ti, observed 
at r = and t — to is given by 

«- ° + (32 ) 



nRi d L 

The angular diameter distance (d A ) is defined as the ratio of the source diameter to its 
angular diameter 

d A = ^ = r 1 R 1 = d L (l + z)- 2 . (33) 



Applying eq.(29) we obtain 



dA = 2^ \TT7] • (34) 
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The maximum of cIa for our model occurs at 

l + ^ = [ 3 (l_ n )]l/(2-3n) (35) 

while for the Freese et al. occurs at [2] 

l + ^=(^) 2/(3 ^ 2) - (36) 
For n = 1/2 we recover the ES result, viz. z m = 5/4. 

7. Number Counts 

The number of astronomical sources in the volume element dV at coordinate r = r 1 is 
given by 

dN = nRhjdttdn (37) 
where n is the number density of the sources and 

dV = R 3 ridndr u dtt = smOdOdcj) (38) 

Since the number of these object is conserved, we can write 

dN = noRlrldttdn (39) 

or 

dN = nv(r l R ) 2 dtt d ^f±dz (40) 
Upon using eq.(29) this transforms into 

For n = 1/2 we recover the ES result, namely 

to. [1 -(! + ,)-*]» 

This has a maximum value at z — [ g~^ ] 1//( - 2 ~ 3n ^ — 1 whereas in the case of ES this occurs 
at z — 16/9 ( i.e. when n = 1/2). 

Equations (17), (23), (29), (34) and (41) are our predictions for the corresponding 
cosmological tests. They are in agreement with the results of the Freese et al. provided 
that the substitution Q = 1 — (3 = 2 — 2nis made. According to the Freese et al. only 
models with /3 < 0.1 are consistent with the Standard Model nucleosynthesis constraints. 
Hence one obtains a restriction on n, viz. n < 0.55 in that epoch. 
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8. Concluding Remarks 

In this work we have discussed the cosmological implications of the bulk viscosity model 
with variable G and A. The results for the cosmological tests are compatible with the 
present observations. The model of the Freese et al. is retrieved from this model for a 
particular choice of n. The Einstein - de Sitter results are obtained for the case n — 1/2. 
While the parameter (3 is a free parameter in the Freese et al. model, it is a measure of 
viscosity in our present model. These tests are found to depend on (3 (or n) only. The 
proposed model resembles the Standard Model in some respect and gives rise to inflation 
for the case n — 1. The model shows the importance of the influence of the bulk viscosity 
in the evolution of the Universe. Further studies are needed to account for the existence 
of this bulk viscosity. 
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